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ABSTRACT 

Context. Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the 

environments of low-mass star forming regions, and in particular the Class protostar IRAS 16293-2422. 

Aims. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of 

the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe 

the fundamental li j - Oo.o transition of the ortho-DiO molecule, inaccessible from the ground, and to determine the ortho-to-para 

D2O ratio. 

Methods. We have detected the fundamental transition of the ortho-DaO molecule at 607.35 GHz towards IRAS 16293-2422. The 

line is seen in absorption with a line opacity of 0.62 ± 0.1 1 (Icr). From the previous ground-based observations of the fundamental 

li,o -lo,i transition of para-D20 seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 ± 0.05 (Icr). 

Results. We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, 

we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 cr level of uncertainty, to be compared with the 

thermal equilibrium value of 2: 1. 

Key words, astrochemistry - ISM: individual (IRAS 16293-2422) - ISM: molecules 



19 
20 
21 



25 



1. Introduction 



2. Observations and results 



Among all molecules in interstellar space, water is a special one 
because of its dominant role in the cooling of warm gas and 
in the oxygen chemistry as well as for its role in the chem- 
istry of the atmospheres of exoplanets and its potential connec- 
tion with life. Water abundance in cold molecular gas is very 
low because it is frozen onto the interstellar grains and forms 
icy mantles around them. Although water can form theoretically 
via gaseous reactions which first form H20^ and H30^ (e.g. 
Rodgers & Charnley 2002), no observational evidence has been 
collected so far It is believed that the major mechanism of wa- 
ter formation is on grain surfaces. One observable that helps to 
discriminate between the various formation mechanisms is the 
abundance of single and double deuterated water with respect to 
the normal isotopologue. Another potential discriminant can be 
the ortho-to-para ratio (OPR), namely the ratio between water 
molecules with different nuclear spins. In fact, since radiative 
and inelastic collisional transitions between the two ortho and 
para states are strongly forbidden, the OPR is set at the moment 
of the water formation and it is changed by nuclear spin reac- 
tions exchange later on. This can occur either in the gas phase by 
reactions with ions in which actual nuclei change places, or on 
the grain surfaces by interaction with electron spins or, perhaps, 
even other nuclear spins (e.g. Le Bourlot et al. 2000, Limbach 
et al. 2006). Although little is known on the spin exchange in 
the gas phase, it is usually assumed that this is a slow process 
and that the OPR is likely to keep memory of the moment of 
its formation. Emprechtinger et al. and Lis et al. in this vol- 
ume report determinations of the water OPR in several environ- 
ments based on new Herschel observations. The doubly deuter- 
ated isotopologue of water, D2O, consists of two species, ortho 
and para with a nuclear spin statistic weight 2:1. So far D2O 
has only been detected towards the solar type protostar IRAS 
16293-2422 (hereafter IRAS 16293), via the observation of the 
fundamental transition of the para-D20 transition at 316.8 GHz 
(see our Figure [T] Butner et al. 2007). The observed line pro- 
file (see Figure |2]i shows a component in emission with a deep 
absorption at the cloud velocity (~ 4 km s"^). The emission 
component has been attributed to heavy water in the hot corino 
of this source where the grain ices are sublimated and released 
into the gas phase (Ceccarelli et al. 2000; Bottinelli et al. 2004), 
based on the detailed analysis of several HDO Unes observed 
in IRAS 16293 (Parise et al. 2005). The absorption component, 
whose linewidth is 0.5 km s ' , is likely due to the foreground gas 
(molecular cloud and cold envelope). Therefore, the absorption 
component provides a straightforward measure of the column 
density of para-D20 in the cold gas surrounding IRAS 16293. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led principal Investigator consortia and with impor- 
tant participation from NASA 



ortho-DgO 



para-DgO 



In the framework of the Key Program CHESS (Ceccarelli et al., 
this volume), we observed the solar type protostar IRAS 16293 
with the HIFI instrument (de Graauw et al., 2010; Roelfsema 
et al., this volume) on board the Herschel Space Observatory 
(Pilbratt et al., 2010). A full spectral coverage of band lb be- 
tween 554.5 and 636.5 GHz was performed on March 2nd 2010, 
using the HIFI Spectral Scan Double Beam Switch (DBS) mode 
with optimization of the continuum. The fundamental ortho- 
D2O (li,i-Oo,o) transition lies in this frequency range, at 607.35 
GHz (see Figure [B- The HIFI Wide Band Spectrometer (WBS) 
was used, providing a spectral resolution of 1.1 MHz (~0.55 km 
s ' at 600 GHz) over an instantaneous bandwidth of 4 x 1 GHz. 
Note that the data are acquired at the Nyquist sampling, there- 
fore, with 0.5 MHz steps. The targeted coordinates were 0-2000 
= 16'' 32'" 22^:75, ^2000 = - 24° 28' 34.2". The beam size at 
610 GHz is about 35", the theoretical main beam (resp. for- 
ward) efficiency is 0.72 (resp. 0.96), and the DBS reference posi- 
tions were situated approximately 3' east and west of the source. 
The data have been processed using the standard HIFI pipeline 
up to level 2 with the ESA-supported package HIPE 3.01 (Ott 
et al. 2010). The 1 GHz chunks are then exported as fits files 
into CLASS/GILDAS formaQfor subsequent data reduction and 
analysis using generic spectral survey tools developed in CLASS 
in our group. When present, spurs were removed in each 1 GHz 
scan and a low order polynomial (< 2) baseline was fitted over 
line-free regions to correct residual bandpass eff'ects. These poly- 
nomials were subtracted and used to determine an accurate con- 
tinuum level by calculating their medians. Sideband deconvo- 
lution is computed with the minimisation algorithm of Comito 
& Schilke (2002) implemented into CLASS using the baseline 
subtracted spectra and assuming side-band gain ratio to be unity 
for all tunings. Both polarisations were averaged to lower the 
noise in the final spectrum. The continuum values obtained are 
well fitted by straight lines over the frequency range of the whole 
band. The single side band continuum derived from the polyno- 
mial fit at the considered frequency was added to the spectra. 
Finally, the deconvolved data were analysed with CASSIS soft- 
warq3 Exact measurements of the main beam efficiency have 
not been performed on planets at the time of our observations. 
However, we are dealing with absorption measurements, and are 
only interested in the relative depth of the absorption compared 
to the continuum level. Consequently we present in the follow- 
ing the spectrum (Figure 2) and parameters (Table 1) in T* for 
the ortho-D20 line. The bottom panel of Figure |2] shows the re- 
sulting HIFI spectrum, with the measured continuum level of 
(234 ±19) mK (where the error includes the statistical error 
only). Note that the absolute calibration doesn't matter for ab- 
sorption, since lines and continuum are affected the same way. 
Therefore, the main source of the uncertainty is the accuracy in 
the continuum. The achieved rms is about 12 mK in T*, in the 
0.5 MHz frequency bin. The fundamental ortho-D20 transition 
at 607349.449 MHz is well detected in absorption against the 
strong continuum, at the velocity of ~ 4 km/s. No other lines in 
the image sideband are expected at this velocity. The parameters 
of the line, obtained using CASSIS, which takes into account the 
ortho and para D2O forms separately from the Cologne Database 
for Molecular Spectroscopy (Miiller et al. 2005, Briinken et al. 
2007), are reported in Table 1 . We report in the same table also 
the parameters of the para-D20 (1 i,o-lo,i) fundamental line pre- 



Fig. 1. Energy levels for the detected fundamental lines of D2O. 
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Table 1. Derived parameters of the ortho and para D2O fundamental lines. Note that the parameters are in T* for ortho-D20 and 
Tmb forpara-D20 (see text). 
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Fig. 2. Profile of the para-D20 (li,o-l(),i) fine (histogram) ob- 
served at JCMT (upper panel), as well as the 3 components gaus- 
sian fit (solid line) and ortho-D20 (1 i,i-lo,i) line observed with 
HIFI (bottom panel). 



viously observed at the JCMT, published by Butner et al. (2007), 
at a rest frequency of 316799.81 MHz. The data were retrieved 
from the JCMT archive and reduced within CLASS. We per- 
formed a 3-component gaussian fit with CASSIS and the result- 
ing fit is reproduced in Figure |2]on top of the data in main beam 
temperatures. The para-D20 line in emission has an intensity 
of 0.10 + 0.02 K in main beam temperature, and a linewidth of 
4.01 ± 0.77 km s^'. The bright line at a V,,, of 10.1 km s"' is 
likely due to CH3OD (see Butner et al. 2007) with an intensity 
of (0.16 + 0.01)K and a Hnewidth of (4.6 + 0.5)km s"'. The pa- 
rameters for the resulting fit of the para-D20 absorption line are 
quoted in Table 1 . 



3. Determination of the D2O OPR 

Crimier et al. (2010) have used the JCMT SCUBA maps of 
IRAS 16293 at 450 yum and 850 //m (and other data) to recon- 
struct the structure of the IRAS 16293 envelope. From this work, 
one can compute the expected continuum in the HIFI beam at 
607 GHz (0-D2O line). Using the SED of Crimier et al. (Fig. 
1 panel d) and their Table 1, the IRAS 16293 flux is 270 + 108 
Jy at 450 fixn and the HIFI beam contains approximately 80% 
of the total source flux (Fig 1, panel b). One can note that the 
SED steep slope yields the flux at 607 GHz to be smaller than 



the one at 450 //m (~ 660 GHz) by about 30%, making the ex- 
pected flux at 607 GHz to be about 0.7 x 0.8 x (270 + 108) Jy 
i.e. (0.34 + 0.14) K, using the HIFI Jy to K conversion factor 
(C. Kramer : Spatial response, contribution to the HIFI frame- 
work document), in perfect agreement with the observed contin- 
uum value (~ 0.33 K in main beam temperature). Most of the 
continuum, more than 70% (resp. 80%) of its peak emission at 
316 GHz (resp. 607 GHz) is emitted from a region of about 900 

AU in radius (~ 15 in diameter). The absorption of the con- 
tinuum by heavy water is most likely due to the cold envelope 
surrounding IRAS 16293 as well as the parent cloud, much more 
extended than the continuum emitting region. Note that, as far 
as the sizes of the absorbing layer are larger than the sizes of 
the region emitting the continuum, the line-to-continuum ratio 
does not depend on the sizes of the telescope beam used for the 
observations. Therefore, we can compute the D2O OPR directly 
from the line-to-continuum ratios of the JCMT and Herschel ob- 
servations, with no further correction. Note also that the para- 
D2O line has an emission component that Butner et al. (2007) 
attributed to the hot corino region, whereas here we are dealing 
with an absorption component only. On the contrary, the ortho- 
D2O line reported here shows an absorption only as the emission 

component is very likely diluted in the 35 HIFI beam, much 

larger than the I5" JCMT beam at 316 GHz. 

Adopting the density and temperature profiles of the enve- 
lope of IRAS 16293 (Crimier et al. 2010), the gas at a distance 
larger (in radius) than 900 AU has a temperature lower than 30 
K and a density lower than about 5 x 10^ cm"^ (see Figure [3). 
Thus, given the temperature of the gas absorbing the D2O lines, 
we only consider the first two levels of each D2O form. We use 
the recently computed collisional rates for the two fundamen- 
tal deexcitation transitions of ortho and para-D20 with para-H2 
in the 10-30 K range: 2.3 x 10"" and 3.8 x 10"" cm^ s"' re- 
spectively (Wiesenfeld, Faure & Scribano, in preparation). At 
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Fig. 3. Density profile (solid line) and gas temperature pro- 
file (dot-dashed) of the IRAS 16293 envelope, as computed by 
Crimier et al. (2010). Values at a radial distance of about 900 
AU are also indicated (see text). A distance of 120 pc has been 
used in this recent computation (Loinard et al. 2008). 
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the low temperatures found in the cold envelope, it is Ukely that 
H2 is mainly in its para form (Pagani et al. 2009, Troscompt et 
al. 2009). With the collisional rates given above, the critical den- 
sities of the ortho- and para- D2O fundamental transitions are 
1x10'^ and 2x10^ cm"^ respectively, and the upper levels of the 
two transitions are only moderately sub-thermally populated for 
a density of 5 x 10^ cm""*. For a two-level system, the species 
column density can be computed as follows: 



N„ 



Sttv^ 



rAV 



V^ Q(Tex) exp(Eu/kTex) 



Auic^ 2V/n2 gu [exp(hv/kTex)- 1] 



(1) 



where A„/ is the Einstein coefficient (2.96 x 10"^ s"' for the ortho 
transition and 6.3 x 10 "* s ' for the para transition), E„ is the up- 
per level energy (E„/k=15.2 K for the para transition and =29.2 
K for the ortho transition), g„ is the upper statistical weight (3 
X (2Jh-1) for the para transition, 6 x (2J+1) for the ortho tran- 
sitions), V is the frequency (316.79981 GHz for the para tran- 
sition and 607.349449 GHz for the ortho transition), AV is the 
linewidth (cm s"') and r is the opacity at the line center T^v is 
the excitation temperature and Q(T(.^) is the corresponding par- 
tition function. In the approximation of the escape probability 
formalism, T^, v is defined by the equation: 



hv/k 



hv/kTk-Hln(l-HAuiy6/Cui) 



(2) 



where Cui - yui x ncoiusion, ncoiiision being the density of the col- 
lision partner (in this case para-H2) and yui being the collisional 
rate in cm-' s"' (values given above). The /3 parameter represents 
the probability that a photon at some position in the cloud es- 
capes the system. For a static, spherically symmetric and homo- 
geneous medium, Osterbrock and Ferland (2006, Appendix 2) 
derives this parameter as a function of the optical depth t in the 
direction of the observer. The opacity at the line center is ex- 
pressed as a function of the line depth (Tabs - Tc-T/,) and the 
continuum (Tc): 



T = -ki 1- 



Tc - Jv(Tex) + Jv(Tcmb)) 



(3) 



Where Jv(Tex) - (hv/k)/(exp(hv/k) - 1) and Tcmb is the cosmic 
microwave background radiation temperature (2.73 K). In the 
limit of T » 1, Tc - Tabs ~ Jv(Tex) - Jv(Tcmb), and Tex ~ 5K. 
Since the D2O transitions are probably optically thin, we can rea- 
sonably assume that Tex is lower than 5 K and Jv(Tex) - Jv(Tcmb) 
is negligible. 

As discussed above, we assume that the absorbing layer is 
much larger than the continuum emitting region. Considering 
the uncertainty on the H2 density (lower than 5 x lO*" cm""*) 
and the kinetic temperature (lower than 30 K), we applied the 
method described above to determine the column densities with 
ne, - lO^cm"-' and Tkin ~ 20K. Table [1] lists the computation 
of the optical depths for both lines as well as the corresponding 
uncertainties. Since t = -ln(TL/Tc) the uncertainty in the line 
optical depth is given by 6t - exp(r) x (5(Tl/Tc). Our compu- 
tation yields an OPR equal to 1.1 ± 0.4 with the corresponding 
column densities: Nonho - (8.7 ± 2.1) 10" cm"^ and Npara - 
(7.8 + 2.6) 10" cm"^. All errors here are 1 cr). Both lines are 
optically thin and their Tp ^ are lower than 5 K. Note that decreas- 
ing the density and/or the kinetic temperature doesn't change the 
OPR by more than 10%. Therefore, the OPR is lower than 2.4 
at a 3 o" level of uncertainty (where we added the 3 cr statistical 
error and the mentioned 10% to the 1.1 value). We assumed (see 
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Fig. 4. Upper limit on the measured D2O OPR (2.6, see text) as 
a grey box and the Boltzmann value (dotted-dashed line) as a 
function of temperature. 



section 2) that the relative gains on the lower and upper side- 
bands are equal. Since we do not have any information concern- 
ing the sideband ratio at the frequency of the D2O line, we can 
only introduce a maximum uncertainty of 16%, corresponding 
to the overall calibration budget for band lb. The resulting up- 
per limit on the OPR is therefore increased to about 2.6. Figure 
IHshows the measured OPR interval and the thermal equilibrium 
as a function of the gas temperature. 



4. Conclusions 

A discussed in |3] the gas absorbing the D2O line lies at more 
than 900 AU from the center and has a temperature lower than 
30 K. The comparison between the upper value of the measured 
D2O OPR and the thermal equilibrium value shows that they are 
consistent with a gas at a temperature larger than about 15 K (at 
a 3 cr level of confidence), and, therefore, with the assumed ab- 
sorbing gas location. On the other hand, the D2O gas could have 
formed in a previous phase, where the gas was colder, and, in this 
case, it means that it had the time to thermalise to the Boltzmann 
value. Unfortunately, given the poor knowledge of the mecha- 
nisms that can exchange the D2O spins (see the Introduction), it 
is difficult to infer here the timescale for this change and, con- 
sequently, to give a lower limit to the object age. On the other 
hand, the relatively large uncertainty in the OPR derived here 
does not allow either to exclude a non- thermal equilibrium sit- 
uation. Higher S/N observations will be needed to lower the un- 
certainty on the OPR value and give a more robust result. 

Using the density and temperature profiles of the envelope of 
IRAS 16293 by Grimier et al. (2010), the column density of the 
gas colder than 30 K is about 1 x 10^^ cm"^. Therefore, the D2O 
abundance (with respect to H2) is about 2x10"". An estimate of 
the water abundance profile will soon be available with the HIFI 
observations with a much higher spatial and spectral resolution 
than the one provided by the ISO observations (Geccarelli et al. 
2000). D2O molecules might form with one OPR, but then could 
freeze out on grains surfaces that could modify the ratio and then 
get desorbed. Due to the high uncertainty in the H2O abundance, 
we cannot at the time being completely exclude or confirm for- 
mation through grain surface chemistry. A modeling of the OPR 
evolution is beyond the scope of the present letter. With an im- 
proved calibration and better understanding of the instrumen- 
tal effects, a more accurate determination of the D2O OPR in 
this source and potentially other sources will be possible. Also, 
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ALMA may hopefully yield an answer in a near future with the 
observation of cold D2O with a higher spatial resolution. 

In summary, this Letter presents the first tentative to estimate 
the OPR for the D2O molecule, demonstrating the outstanding 
capabilities of the HIFI instrument. The poor knowledge of the 
mechanisms of exchange of the nuclear spins and the relatively 
large error in the derived OPR prevent to drawing firm conclu- 
sions on the formation of heavy water at that time. 
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